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ABSTRACT 

Context. Stellar shells observed in many giant elliptical and lenticular as well as a few spiral and dwarf galaxies presumably result 
from galaxy mergers. Line-of-sight velocity distributions of the shells could, in principle, if measured with a sufficiently high signal- 
to-noise ratio, constitute a method to constrain the gravitational potential of the host galaxy. 

Aims. Merrifield & Kuijken (1998) predicted a double-peaked line profile for stationary shells resulting from a nearly radial minor 
merger. In this paper, we aim at extending their analysis to a more realistic case of expanding shells, inherent to the merging process, 
whereas we assume the same type of merger and the same orbital geometry. 

Methods. We used an analytical approach as well as test particle simulations to predict the line-of-sight velocity profile across the 
shell structure. Simulated line profiles were convolved with spectral PSFs to estimate peak detectability. 

Results. The resulting line-of-sight velocity distributions are more complex than previously predicted due to nonzero phase velocity 
of the shells. In principle, each of the Merrifield & Kuijken (1998) peaks splits into two, giving a quadruple-peaked line profile, which 
allows more precise determination of the potential of the host galaxy and contains additional information. We find simple analytical 
expressions that connect the positions of the four peaks of the line profile and the mass distribution of the galaxy, namely, the circular 
velocity at the given shell radius and the propagation velocity of the shell. The analytical expressions were applied to a test-particle 
simulation of a radial minor merger, and the potential of the simulated host galaxy was successfully recovered. Shell kinematics can 
thus become an independent tool to determine the content and distribution of the dark matter in shell galaxies up to ~100kpc from 
the center of the host galaxy. 

Key words. Galaxies: kinematics and dynamics - Galaxies: halos - Methods: analytical - Methods: numerical - Galaxies: elliptical 
and lenticular, cD - Galaxies: interactions 



1. Introduction 

Several methods have been used to measure the gravitational 
potentials and their gradients of elliptical galaxies, including 
strong and weak gravitational lensing (e.g., Gavazzi et al.| 2007; 
Mandelba um et al.|2008[|Auger et al.|2010 1, X-ray observations 
of hot gas in the massive gas-rich galaxies (e.g.,[Fukazawa et al ' 



2006, Na gino & Matsus hita 2009[ |Churazov et al.| |2008 Das 
et al.|2010| i, rotation curves of detected disks and rings of neutral 



hydrogen (e.g., |Weijmans et al.||2008| l, stellar-dynamical mod 



eling from integrated light spectra (e.g., Weijmans et al.| 2009 



de Lorenzi eTaT1[2009l |Churazov et al.||20101 [Thomas et al.| 
201 IK and the use of tracers such as planetary nebulae, glob- 
ular clusters, and satellite galaxies (e.g., |Coccato et al.||2009 



Nierenberg et al.|2011||Deason et al.|2012||Norris et al.|20T2j . 
All these methods have their limits, such as the redshift of the 
observed object, the luminosity profile, and gas content. In par- 
ticular, the use of stellar dynamical modeling is plausible in the 
wide range of galactic masses, as long as spectroscopic data are 
available. However, it becomes more challenging beyond a few 
optical half-light radii. Other complementary gravitational trac- 
ers or techniques are required to derive mass profiles in outer 



parts of the galaxies. When comparing independent techniques 
for the same objects at similar galactocentric radii, discrepancies 
in the estimated circular velocity curves were revealed together 
with several interpretations (e.g., |Churazov etaL[2 010; Das et al. 
|2010[ ). The compared techniques usually employ modeling of 
the X-ray emission of the hot gas (assuming hydrostatic equilib- 
rium) and dynamical modeling of the optical data in the massive 
early-type galaxies. Therefore, even for the most massive galax- 
ies with X-ray observations available, there is a need for other 
methods to independently constrain the gravitational potential at 
various radii. 

Shell galaxies are galaxies that contain arc-like fine features, 
which were first noticed by Arp [ ( |1966| l. These structures are 
made of stars and form open, almost concentric arcs that do not 
cross each other. Shells are relatively common in elliptical or 
lenticular galaxies. At least 10% of all these galaxies in the lo- 
cal universe possess shells. Nevertheless, shells occur markedly 
most often in regions of low galaxy density, and perhaps up to 
half of E and SO galaxies in these environments are shell galax- 



ies (Malin & Carter 


1983 


|Schweizer||1983||Schweizer & Ford| 


1985, Colbert et al. 


2001 


i. Shells can also be associated with 
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dust ( [Sikkema et al.||2007| [Stickel et al.||2004|) and neutral hy- 



drogen emission ( Schiminovich et al. 
ISancisi 1996; Petric et a: 



tion, Charmandaris et al. 



molecular gas in the shel 



1994 1995 Balcells & 



latter case also X-ray ones) in a few shell galaxies. The results 
support accretion events in the history of shell galaxies. 



1997||Horeiiou et al.|2001) . In addi- Merrifield & Kuijken (1998), hereafter MK98 studied the 



(2000) detected the presence of dense 



sofNGC5128. 



Shells are thought to be by-products of minor mergers of 
galaxies (Quin n|1984|), although they c an also be formed during 
major mergers ( |Her nquist & Spergel 1992). The most regular 
shell systems are believed to result from nearly radial mergers 
( [Dupraz & Combes|[T986l |Hernquist & Quinn|[T988] ). When a 
small galaxy enters the sphere of influence of a big elliptical 
galaxy on a radial or close-to-radial trajectory, it disintegrates 
and its stars begin to oscillate in the potential of the big galaxy. 
At their turning points, the stars have the lowest speed and thus 
tend to spend most of the time there, where they pile up and 
produce arc-like structures in the luminosity profile of the host 
galaxy when viewed perpendicular to the axis of the collision. 

Measurement of the number and distribution of shells can, in 
principle, yield to an approximate estimate of the mass distribu- 
tion of the host galaxy and the time since the merg er ([Q uinn 



1984 



Cana 



Dupraz & Combes|[T9"8"6} |Hernquist & Quinn| [T987a b 



izo et al. 1 12007] ). But both of these observables are sen 



sitive to details such as the dynamical friction and the gradual 
decay of the cannibalized galaxy during the merger (Dupraz & 



Combes 1987; James & Wilkinso n]T987l|Heisler & White| 1990 
Ebrova et al. 11201 0b. Moreover, if the core of the cannibalized 



galaxy survives the merger, new generations of shells are added 
during each successive passage. This was predicted by Dupraz 
& Combes ( |1987| ) and successfully reproduced by Bartoskova 
et al.| ( |20lT |) in self-consistent simulations. All these effects com- 



plicate the simulations to such an extent that the interest in shell 
galaxies largely faded by the end of the 1980s. Recently, this 
topic has raised interest again, thanks to the discovery of shells 
in a quasar host galaxy (Canalizo et al. 2007) and shell struc- 
tures in M31 (Fardal et al. 2007 2008) and in the Fornax dwarf 



( Coleman et al.|2004[ ). |Helmi et al.| ( |2003[ ) suggested that ring- 
like stellar structures, including the one observed in the outer 
disk of the Milky Way (the so-called Monoceros ring), could be 
analogous to shells. A significant number of shells is also con- 
tained in the early-type galaxy sample of the ongoing ATLAS 3D 
project, including images of galaxi es with a surface brightness 
down to 29mag/arcsec 2 



Krajnovic et al. 2011 



Due 



(see, e.g. 

|et al.||201 lb . |Kim et al.| ( pOT2] ) identified shells in about 6% of 
a sample of 65 early-type galaxies from the Spitzer Survey of 
Stellar Structure in Galaxies (S 4 G). Shells also appear to be suit- 
able for indirect detection of dark matter via gamma-ray emis- 
sion from dark matter self-annihilations (Sanders on et al.|2012| ). 
About 70% of a complete sample of nearby (15-50 Mpc) lumi- 
nous (Mb < -20 mag) elliptical galaxies were found to show 
tidal features by Tal et al. (2009). Faint structures, including 
shells and other signatures of recent gravitational interaction 
(tidal tails and streams), were found in the Sloan Digital Sky 
Survey (SPSS). [KavTrajl pOlO] ) identified 18% of early-type 
galaxies (ETGs) in the SDSS Stripe82 sample as having dis- 
turbed morphologies; similarly, Miskolczi et al.| ( |201 1) found 
tidal features in 19 % of their sample of galaxies from SDSS 
DR7. Observations of warm gas by Rampazzo et al. (2003 1 in 
five shell galaxies showed irregular gaseous velocity fields (e.g., 
a double nucleus or elongated gas distribution with asymmetric 
structure relative to the stellar body), and in most cases, gas and 
stellar kinematics appe ar decoupled. |Rampazzo et al.| ( |2007| ), 



Marino et al. ( |2009[ ), and |Trinchieri et al . (2008 ) investigated star 
formation histories and hot gas content using the NUV and FUV 
Galaxy Evolution Explorer (GALEX) observations (and in the 



oretically the kinematics of a stationary shell, a monoenergetic 
spherically symmetric system of stars oscillating on radial orbits 
in a spherically symmetric potential. They predicted that spec- 
tral line profiles of such a system exhibit two clear maxima, 
which provide a direct measure of the gradient of the gravita- 
tional potential at the shell radius. The first attempt to analyze 
the kinematical imprint of a shell observationally was made by 
|Romanowsky et aT1 ( 2012| ), who used globular clusters as shell 
tracers in the early-type galaxy M87. Far dal et al.| ( |2012| ob- 
tained radial velocities of giant stars in the so-called western 
shelf in M3 1 Andromeda galaxy. They successfully analyzed the 
shell pattern in the space of velocity versus radius. 

Nevertheless, real-world shells are not stationary features. 
The stars of the satellite galaxy have a continuous energy distri- 
bution, and so at different times, the shell edge is made of stars 
of different energies, as they continue to arrive at their respec- 
tive turning points. Thus, the shell front appears to be traveling 
outwards from the center of the host galaxy and shell spectral- 
line profiles ar e more complex (|Jflkova et al. 2010 Ebrova et al. 
|20TT] see also |Fardal et al.|2012| ). 

In this paper, we derive spectral-line profiles of nonstation- 
ary shells. We assume that shells originate from radial minor 
mergers of galaxies, as proposed by |Quinn (1984]). We find that 
both of the original MK98 peaks in the spectral line are split 
into two, resulting in a quadruple shape, which can still be used 
to constrain the host galaxy potential and even bring additional 
information. We outline the simplified theoretical model and de- 
rive the shell velocities in Sect. |2) and describe the origin of 
the quadruple line profile in Sect.]3| In Sect. [4] we derive equa- 
tions connecting the observable features of the quadruple-peaked 
line-of-sight velocity distribution (LOSVD) with parameters of 
the host galaxy potential in the vicinity of the shell edge. We 
compare these analytical predictions with the theoretical model 
(Sect. [5]) and with results of test-particle simulations of the radial 
minor merger (Sect.|6|. Section [6] also demonstrates the deriva- 
tion of the galactic potential from the simulated spectral data. 



2. Model of radial oscillations 

If we approximate the shell system with a simplified model, we 
can describe its evolution completely depending only on the po- 
tential of the host galaxy. The approximation lies in the numeri- 
cal integration of radial trajectories of stars in a spherically sym- 
metric potential. The distribution of energies of stars is contin- 
uous, and these stars were released from a small volume in the 
phase space. We call it the model of radial oscillations, and it cor- 
responds to the notion that the cannibalized galaxy came along a 
radial path and disintegrated in the center of the host galaxy. As 
a result the stars were released at one moment in the center and 
began to oscillate freely on radial orbits. This approach was first 
used by |Quinn|(|1984[), followed by [Dupraz & Combes| (1986 
T987] ) and pernquist & Qurnn[ ( 11987a|b] r~ 



2. 1 . Turning point positions and their velocities 

In shell galaxies, the shells are traditionally numbered accord- 
ing to the serial number of the shell, n, from the outermost to 
the innermost (which in the model of radial oscillations for a 
single-generation shell system corresponds to the oldest and the 
youngest shell, respectively). If the cannibalized galaxy comes 
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from the right side of the host galaxy, stars are released in the 
center of the host galaxy. After that, they reach their apocenters 
for the first time. But a shell does not form here yet, because the 
stars are not sufficiently phase wrapped. We call this the zeroth 
oscillation (the zeroth turning point) as we try to match the num- 
ber of oscillations with the customary numbering scheme of the 
shells. We label the first shell that occurs on the right side (the 
same side from which the cannibalized galaxy approached) with 
n = 1. Shell no. 2 appears on the left side of the host galaxy, 
no. 3 on the right, and so forth. 

In the model of radial oscillations, the shells occur close to 
the radii where the stars are located in their apocenters at a given 
moment (the current turning point, rrp, in our notation). The 
shell number n corresponds to the number of oscillations that 
the stars near the shell have completed or are about to complete. 
The current turning point rrp must follow the equation 



f = (« + l/2)r(r TP ), 



(1) 



where f is the time elapsed since stars were released in the cen- 
ter of the host galaxy. T(r) is the period of radial motion at a 
galactocentric radius r in the host galaxy potential </>(r): 



7» = V2 f [4>(r)-cf>(r')Y l/2 dr'. 
Jo 



The position of the current turning point evolves in time with 
a velocity given by the derivative of Eq. ([T]i with respect to radius 



v TP (r; n) = dr/df : 



1 



1 



dt/dr n + 1/2 



(dT(r)/drT 



(3) 



We can clearly see from this relation, which was first derived by 
Quinn ( 1984 1, that any further turning point (turning point with 
higher n) at the same radius moves more slowly than the former 
one. Thus causes a gradual densification of the space distribution 
of the shell system with time. 

Technically, the reason for this densification is that the time 
difference between the moments when two stars with similar en- 
ergy reach their turning points is cumulative. Let Af be the dif- 
ference in periods at two different radii r a and r\, (with r a < r\,, 
on the right). The radius where stars complete the first oscilla- 
tion moves from r a to r\, in Af. But in the second orbit on the 
left, the stars from r\, will already have a lag of Af behind those 
from r a and will just be getting a second one, so the third one 
(the second on the same side) reaches from r a in 2 x Af. Every 
nth completed oscillation on the right side, then moves n times 
more slowly than the first one. The situation is similar on the left 
side, and the shell system is getting denser. Moreover, the turn- 
ing point has an additional lag of 1 /2T(rjp), because the stars 
were released in the center of the host galaxy before their zeroth 
oscillation. This is the source of the factor (n + 1/2) in Eqs. ([T]i 
and @. 



to the plane of the sky (Nulsen 1989). For the nth shell, this is 
the largest radius where stars about to complete their nth oscilla- 
tion are currently located. This radius is always larger than that 
of the current turning point of the stars that are completing their 
nth oscillation. Thus, the shell edge consists of outward-moving 
stars about to complete their nth oscillation. 

|Dupraz & Combes| fl986) state that the stars forming the 
shell move with the phase velocity of the shell. While we show 
that this holds only approximately, we use this equality in Sect. [4] 
to derive the relation between the shell kinematics and the poten- 
tial of the host galaxy. 

The position of a star, r„, at a given time f since the release 
of the star in the center of the host galaxy is given by an implicit 
equation for r* and is a function of the star energy, or equiva- 
lently the position of its apocenter r ac .QFor stars with the integer 
part of f/[27\r ac )] odd, the equation reads: 



f = (n + l)V2 £" [^(r ac )-^(r')r 1/2 dr'- 
- /;* [2Wr ac )-<Xr'))r 1/2 dr'. 



(4) 



For stars that have completed an even number of half -periods 
(only such stars are found on the shell edge), the equation is 



(2) f = „V2 £' [ ( p(r ic )- ( p(r')]- 1 ' 2 dr'+ 



Jo 
+ Jo 







[2(<Kr ac )-0(r'))r 1/2 dr'. 



(5) 



The first term in Eq. |5]) corresponds to n radial periods for the 
star's energy (n is maximal so that nr(r ac ) < f), while the other 
term corresponds to the time that it takes to reach radius r* from 
the center of the galaxy. Even for the simplest galactic potentials, 
these equations are not analytically solvable and must be solved 
numerically. 

The position of the nth shell r s equals the maximal radius 
r*,max that solves Eq. Q for the given n. The shell velocity v s is 
obtained from the numerical derivative of a set of values of r, max 
for several close values of f . 

The stellar velocity at the shell edge is obtained by inserting 
7"*,max with its corresponding r ac into: 



(6) 



v(r*) = + ^2[(p(r ac ) - 0(r„)]. 

For the stars following Eq. ([5]), the velocity will be positive; for 
the rest, it will be negative. 

It is clear that v(r* ;lnax ) < v s . Actually, v(r tiinax ) is always 
slightly lower than the phase velocity of the shell (Table [TJ. 
Meanwhile, the position of the shell for a given time is not 
far from the current turning point, and their separation changes 
slowly. Thus, the velocity of the turning points given in Eq. Q is 
a good approximation for the shell velocity (Fig. [3}. Equation ([3]l 
is not generally solvable analytically either, but the numerical 
calculation of vtp is much easier than determining the true ve- 
locity v s as described in this chapter. 



2.2. Real shell positions and velocities 

Even in the framework of the radial oscillation model, the posi- 
tion and velocity of the true edge of the shell cannot be expressed 
in a straightforward manner. Photometrically, shells appear as a 
brightening in the luminosity profile of the galaxy with a sharp 
cut-off. This is because the stars of the cannibalized galaxy oc- 
cupy a limited volume in the phase space. With time, the shape 
of this volume gets thinner, more elongated, and wrapped around 
invariant surfaces defined by the trajectories of the particles, 
increasing its coincidence with these surfaces. A shell appears 
close to the points where the invariant surface is perpendicular 



2.3. Kinematics of shell stars 

In the same model, we can also describe the LOSVD of a shell 
at a given time f, for a given potential of the host galaxy </>(r). In 
this paper, we model the host galaxy potential of the host galaxy 
as a double Plummer sphere, as described in Sect. 



6.1 



Eqs. Q and |5]l give the actual star position r* and the shell 
number n for any apocenter r ac in a range of energies. The radial 



1 We denote the apocenter of the star corresponding to its energy as 
r ac , whereas r TP (the current turning point) is the radius at which the 
stars reach their apocenters at the time of measurement. 
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far half of the galaxy 
near half of the galaxy 
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A+ 



observer 





j [km/s] 



MK98 



for a stationary shell. Left: Scheme of the kinematics of a shell 



Fig. 1. Kinematics of a moving shell. Compare with Fig. 3 in ! 

with radius r s and phase velocity v s . The shell is composed of stars on radial orbits with radial velocity v r and LOS velocity vi 
Right: The LOSVD at projected radius R = 0.9r s , where r s = 120 kpc (parameters of the shell are highlighted in bold in Table[T]l, 
in the framework of the model of radial oscillations (Sect. [23). The profile does not include stars of the host galaxy, which are not 



part of the shell system, and is normalized, so that the total flux equals one. (a) The LOSVD showing separate contributions from 
inward and outward stars; (b) the same profile, separated for contributions from the half of the host galaxy closer to the observer 
(the one including point B) and the more distant half (includes point A). 



velocity of a star on the particular radius is given by inserting 
the corresponding pair of r ac a r, in Eq. (|6j. Naturally, the pro- 
jections of these velocities to the selected line-of-sight (LOS) 
form the LOSVD. To reconstruct the LOSVD, we have to add 
an assumption about the behavior of shell brightness in time. In 
Sect. |4.2| we show our choice of the behavior and also illustrate 
that the particular choice does not matter much. 



3. Quadruple-peaked LOSVD 

Fig. [T] illustrates a measurement of the LOSVD of stars in the 
shell, which is composed of inward and outward stars on radial 
trajectories. The stars near the edge of the shell move slowly. 
But it is clear from the geometry that contributions add up from 
different galactocentric distances, where the stars are either still 
traveling outwards to reach the shell or returning from their 
apocenters to form a nontrivial LOSVD. MK98 showed that the 
maximal contribution to the LOSVD comes from stars at two 
particular locations along the line of sight (A and B), both of 
which are at the same galactocentric distance. 

In |MK 98 s stationary shell model, inward stars at the same 
radius differ from outward stars only in the sign of the LOS ve- 
locity vios. This is not true when the edge of the shell moves out- 
wards with velocity v s . At any given instant, the stars that move 
inwards are returning from a point where the shell edge was at 
some earlier time, and so their apocenter is inside the current 
shell radius r s . Similarly, the stars that move outwards will reach 
the shell edge in the future. Consequently, the stars that move 
inwards are always closer to their apocenter than those moving 
outwards at the same radius, and their velocity is thus smaller. 



The inward stars move toward the observers in the farther of 
the two MK98] points (A) and away from them in the nearer 
point (B), while the stars moving outwards behave in the op- 
posite manner. Together, there are four possible velocities with 
the maximal contribution to the LOSVD, resulting in its sym- 
metrical quadruple shape shown in Fig.[T] In fact, for a moving 
shell, points A and B are not at the same galactocentric radius for 
inward and outward stars. For inward stars, points A and B are a 
little closer to the center as indicated in Fig.[T] This is discussed 
in Sect.^" 



In the right-hand panel of Fi g.fT] we used the model of radial 
oscillations as described in Sect.|2.3|to illustrate individual con- 



tributions to the LOSVD. MK98| constructed an analytical func- 
tion describing the LOSVD close to the edge of their stationary 
shell model. This function exhibits intensity maxima that coin- 
cide with maximal/minimal velocity, leading to the symmetrical 
double-peaked profile with peaks at the edges of the LOSVD. 
This cannot be shown for a general moving shell, but Fig. [Tj 
demonstrates that the intensity maxima coincide with velocity 
extremes for separate contributions to the LOSVD. 

The separation in velocity between peaks for a given pro- 
jected radius R is given by the distance of R from the edge of 
the shell r s . The profile shown in Fig.[T]corresponds to projected 
radius R = 0.9r s . The closer to the shell edge, the narrower the 
profile is. The separation of the peaks at a given R depends on 
the phase velocity of the specific shell, near which we observe 
the LOSVD. This velocity is, for a fixed potential, given by the 
shell radius and its serial number (Sect. |2.1| i. These effects are 
illustrated in Fig. [2] where we show the positions of the LOSVD 
peaks for the first shell at different radii r s and for a shell at 
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Table 1. Parameters of shells for which the LOSVD intensity 
maxima are shown in Fig. [2] 



Fig. 2. Locations of peaks of the LOSVDs in the framework of 
the model of radial oscillations (Sect. [23) : (a) for the first shell 
at different radii, (b) for the first to the fourth shell at the radius 
of 120 kpc. Parameters of all shells are shown in Table [JJ For 
parameters of the host galaxy potential, see Sect. 6. 1 



120 kpc with different serial numbers n. Note that the higher the 
serial number n at a given radius, the smaller is the difference in 
the phase velocity between the two shells with consecutive se- 
rial numbers and thus in the positions of the respective peaks. 
Parameters of the corresponding shells can be found in Table [JJ 
The radial dependence of the phase velocity of the first four 
shells in the whole host galaxy is shown in Fig. [3] Using Eq. Q, 
we see that the velocity of each subsequent shell differs from the 
first one only by a factor of 3/(1 + 2n). The large interval of the 
galactocentric radii where the shell velocity increases is caused 
by the presence of the halo with a large scaling parameter. In 
fact, we do not show shell velocity, but the velocity of the turn- 
ing points at the same radius. Nevertheless, these are sufficiently 
close. Black crosses show the true velocity of the first shell cal- 
culated for several radii according to the method described in 



Sect. 2.2 For shells of higher n, these differences between the 
phase velocity of a shell and the corresponding turning point 
with consecutive serial numbers are even smaller. 

The radius of a stationary shell is the same as the radius of 
the apocenter of stars (as they all have the same energy), while 
the edge of a moving shell is at the radius, which is always 
slightly further from the center than the current turning points. 
This difference creates an intricate zone between the radius of 
the current turning points and the radius of the edge, where all 
the stars of a given shell move outwards. When the LOS ra- 
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Notes, t: time since the release of stars at the center of the host galaxy, 
in which the shell has reached its current radius calculated in the frame- 
work of the model of radial oscillations (Sect. [2); n: serial number of 
shell (Sect. [2~T| >; r s : shell radius; v s : shell phase velocity according to 
the method described in Sect. |2.2| r TP : galactocentric radius of current 
turning points of the stars at this time given by Eq. v(r»_ mllx ): radial 
velocity of stars at the shell edge; v TP : phase velocity of current turning 
point according Eq. l[3j; v c : circular velocity at the shell edge radius. 
For parameters of the host galaxy, see Sect. |6. 1| The shell that is used 
in Figs.|4]|6]|7] and[5]is highlighted in bold. 




Fig. 3. Dependence of the phase velocity of the turning points 
on the galactocentric radius for the first four shells according to 
Eq. ([3]). For parameters of the host galaxy potential, see Sect. 6. 1 



Black crosses show the true velocity of the first shell calculatec 



for several radii according to the method described in Sect. 2.2 
In fact, the turning point responsible for the current location of 
shell is not at the same radius as the shell edge at the same time, 
but the difference is small (Table [JJ). 

dius from lower radii gets near to the turning points of the stars, 
the inner maxima of the LOSVD approach each other until they 
merge and finally disappear (Fig.Q. We actually see a minimum 
in the middle of the LOSVD closer to the shell edge than the cur- 
rent turning points. The intricate zone is much larger for the first 
shell. For the shell radius of 120 kpc in our host galaxy potential, 
it occupies 6 kpc for the first shell, 2 kpc for the second one, and 
less than one kpc for the fourth shell (Table [TJ. 

4. Relating observables to circular and shell 
velocities 

The nonzero velocity of the shell complicates the kinematics of 
shells in two aspects mentioned above. Due to the energy differ- 
ence between inward and outward particles at the same radius, 
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Fig. 4. Evolution of the LOSVD near the shell edge for the sec- 
ond shell at r s = 120 kpc (parameters of the shell are high- 
lighted in bold in Table [T} for the projected radius 116, 117, 
118, and 1 19 kpc in the framework of the model of radial oscil- 
lations (Sect. 2.3 1. In this model, the current turning points of 
the shell particles are at rtp = 1 17.9 kpc. Beyond this radius, the 
inner maxima disappear. Profiles do not include stars of the host 
galaxy, which are not part of the shell system and are normal- 
ized so that the total flux equals one. For parameters of the host 
galaxy potential, see Sect. RTT 



the LOSVD peak is split into two and the shell edge is not at the 
radius of the current turning point, but slightly further from the 
center of the host galaxy. In this section, we describe the LOSVD 
of such a shell in the approximation of a locally constant galac- 
tic acceleration and shell velocity. In addition, we assume that 
the velocity of stars at the edge of the shell is equal to the phase 
velocity of the shell. 

4.1. Motion of a star in a shell system 

The galactocentric radius of the shell edge is a function of time, 
r s (t), where t — is the moment of measurement and r s (0) = r s o 
is the position of the shell edge at this time. We assume that the 
stars are on strictly radial orbits, that there is a locally constant 
value of the radial acceleration oq in the host galaxy potential 
and a locally constant velocity of the shell edge v s , and that the 
stars at the shell edge have the same velocity as the shell. The 
galactocentric radius of each star is at any time r{t), while f s is 
the time when the star could be found at the shell edge r s (f s ). 
Then the equation of motion and the initial conditions for the 
star near a given shell radius are 



d 2 r(t) 
dt 2 



dr(t) 



dt 

r(f s ) = r s (f s ) = v s t s + r s0 . 
The solution of these equations is 
r(f) = ao(t - f s ) 2 /2 + v s (t - t s ) + r s (f s ), 
v(f) = v s + a (t - f s ), 



(7) 



(8) 



(9) 



(10) 



and the actual position of the star r(0) and its radial velocity v(0) 
at time of measuring (t = 0) are 



r(0) = t 2 s a /2 + r s 



s<)> 



v(0) = v s -a Q t s . 

Eliminating f s from the two previous equations, we get 
v(0) ± = v s + v c V2(l-r(0)/r s0 ), 



(12) 
(13) 



(14) 



where v c = y-ao>" s o is the circular velocity at the shell edge 
radius. 



4.2. Approximative LOSVD 



The projection of the velocity given by Eq. ( 14 1 to the LOS at a 
projected radius R will be 



Vlo 



= VT-* 2 /(r(0)) 2 v(0) ± = 

= Vl - R 2 / (r (0)) 2 [v. ± v c V2 (1 - r(0)/r s0 )] . 



(15) 



Using this expression, we can model the LOSVD at a given pro- 
jected radius for a given shell. For the proper choice of a pair of 
values v c and v s , we can find a match with observed and modeled 
peaks of the LOSVD. 

To model the LOSVD in both frameworks, the model of ra- 



dial oscillations (Sect. |2.3[ ) and the approximative LOSVD by 
Eq 
the s 



( p3| ), we have to add an assumption about the behavior of 
lell brightness in time or in space (as the shell expands with 
time). This behavior depends on the parameters of the merger 
that has produced the shells. It is determined by the energy dis- 
tribution of stars of the cannibalized galaxy in the instant of its 
decay in the center of the host galaxy. For simplicity, we choose 
the density at the surface of a sphere of shell edge radius r s 
to be E sp h(r s (f)) ~ l/r 2 (f), corresponding to a shell containing 
the same number of stars at each moment. The relation between 
S sp h(r s (f)) and the projected surface density near the shell edge 

on the sky Zi os (r s (f)) is Li os (r s (r)) ~ y/r s (t)I. svb (r s (t)). It turns out 
that no reasonable choice of this function has an effect on the 
general characteristics of the LOSVD and the principles of for- 
mation that we describe in this paper. For illustration, we demon- 
strate the LOSVD of S sp h increasing as r 2 and E sp h decreasing as 
1 /r 2 in Fig. [5] For the profiles shown, the ratio of the inner and 
outer peaks changes with the change of the £ sp h, but the peak 
positions are unaffected and the overall shape of the profile does 
not alter significantly. For shells that were created in a radial mi- 
nor merger, we can expect a sharp rise in shell brightness near 
the center of the host galaxy, followed by an extensive area of its 
decrease. The fact that the main features of the LOSVD do not 
depend on the choice of E sp h means that our method of measur- 
ing the potential of shell galaxies is not sensitive to the details of 
the decay of the cannibalized galaxy. 

4.3. Radius of maximal LOS velocity 

MK98 proved that near the edge of a stationary shell, r s , the 



maximum intensity of the LOSVD is at the point where the max- 
imal absolute value of the LOS velocity is. They also proved that 
the maximal absolute value of the LOS velocity vi os max comes 
from stars at the galactocentric radius 



rvmax = ^(R + r s o), 



(16) 



(11) at each projected radius R. 
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Fig. 5. LOSVD of the second shell at r s = 120kpc (parameters 
of the shell are highlighted in bold in Table [TJ for the projected 
radius 108 kpc in the framework of the model of radial oscilla- 
tions (Sect. |2.3) , where the density at the surface of a sphere of 
shell edge radius r s is £ sp h(r s (f)) ~ r 2 (t) for the blue curve and 

2sph('"s(0) ~ l/ r s (0 f° r tne re d one - The profile does not include 
stars of the host galaxy, which are not part of the shell system 
and are normalized, so that the total flux equals one. 
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Fig. 6. Galactocentric radii r lmax that contribute to the LOSVD 



the maximal velocities according to Eq. (16i, which was used 
in the derivation of the approximative maximal/minimal LOS 
velocities (Sect. [5] point [3]) - orange curve, according to the ap- 
proximative LOSVD (Sect. [5] point [2]i - purple curves, and ac- 
cording to the model of radial oscillations (Sect. [23] )- light blue 
curves for the second shell at 120 kpc (parameters of the shell 
are highlighted in bold in Table [T}. For parameters of the host 
galaxy potential, see Sect. |6.1| 



For a moving shell, analogous equations are significantly 
more complex and a similar relation cannot be easily proven. 
Nevertheless, when we apply both results of MK98 we can show 



+ n 



os.max- 



in examples (Figs. [7] [81 |TT] and 12 1 that their use is valid, even 
for nonstationary shells. In the framework of the radial oscil- 
lations model (Sect. |2.3| l, we have shown that the peaks of the 
LOSVD occur fairly close to the edges of distributions of inward 
and outward stars (Fig. |TJ. The peaks are also near the edges of 
the LOSVD, if we divide the LOSVD into the contributions of 
the near and the far half of the galaxy as in Fig.[T](b). The inner 
peak corresponds to inward-moving stars and the outer one to 
outward-moving ones. This approach is used in the equations in 
Sect. 4.4 The maximal LOS velocity corresponds to the outer 
peak and the minimal to the inner one. Reasons and justification 
for use of Eq. ( 16 1 for are discussed in Sect. [5] point|3](see 
also Fig. |6j. 

4.4. Approximative maximal LOS velocity 



Using the results of MK98 we derive an expression for the 
maxima/minima of the LOS velocity corresponding to locations 
of the LOSVD peaks in observable quantities (i.e., the max- 
ima/minima of the LOS velocity, the projected radius, and the 
shell radius) by substituting r vm . dx given by Eq. ( 16 1 for r(0) in 
Eq.fSJ 



^los,max± 



(v s ± I 



VI -fl/rso) 



•4CR/r s o) 2 (l+/?/r s0 )- 2 . 



For the measured locations of the LOSVD peaks Vi os ,max+, 
vios.max-, projected radius R, and shell edge radius r s o, we can 
express the circular velocity v c at the shell edge radius and the 
current shell velocity v s by using inverse equations: 



plos,max+ ^'los,max- 



2^/(1- R/r s0 ) [ 1 - 4 (R/r s0 ) 2 (1 + R/r s0 T 2 ] 



2Vl-4(/?/r s „) 2 (l+/?/r s0 r 2 



(19) 



Alternatively, the value of the circular velocity v c at the shell 
edge radius could be inferred from measurements of positions of 
peaks at two or more different projected radii for the same shell: 
let Avi os = Vi os ,max+ - iVs.max-, where Vi os ,max± satisfy Eq. ([17 I. 
Then, in the vicinity of the shell edge, 



Avios = 2v c ^(R/r s0 - l)[l -4(/?/r s0 ) 2 (l +R/r s0 T 2 ] 
* 2(1 -R/r s0 )v c , 



(20) 



and taking the derivative with respect to the projected radius 



dAv 



los 



dR 



(21) 



which happens to be the same expression as equation (7) in 
MK98 Nevertheless, in MK98] Avi os is the distance between the 
two LOSVD intensity maxima of a stationary shell, whereas in 
our framework, it is the distance between the outer peak for pos- 
itive velocities and the inner peak for negative velocities or vice 
versa. This equation allows us to measure the circular velocity in 
shell galaxies using the slope of the LOSVD intensity maxima 
in the R X vi os diagram. 



(17) 5. Comparison of models 



In this section, we compare three different approaches to the 
theoretical calculation of the maximal/minimal LOS velocities, 
which are equivalent to the positions of LOSVD peaks: 



1. Using the model of radial oscillations as described in 
Sect. |2.3| (these results are plotted with light blue curves in 
relevant figures). This model requires thorough knowledge 
(18) of the potential of the host galaxy, obviously unavailable for 
real galaxies. 
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Fig. 7. LOSVD peak locations for the second shell at the radius 
of 120 kpc (parameters of the shell are highlighted in bold in 
Table [T} according to the approximative maximal LOS veloc- 
ities (Sect. [5] point [3]) given by Eq. (17i (orange curves); the 
approximative LOSVD (Sect. [5] point|2j) given by Eq. (fl~5| > (pur- 
ple curves); and the model of radial oscillations (Sect. [23} (light 
blue curves almost merged with the purple ones). The red line 
shows the position of the LOSVD from Fig. [HJ the black one 
shows the position of the current turning points. For parameters 
of the host galaxy potential, see Sect. 6.1 
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Fig. 8. LOSVD of the second shell at r s = 120 kpc (parame- 
ters of the shell are highlighted in bold in Table [TJ for the pro- 
jected radius R = 0.9r s = 108 kpc according to the approxima- 
tive LOSVD (Sect. [5] point [2jl given by E q. (fT5) > (purple curve) 
and the model of radial oscillations (Sect. 2.3 1 (light blue curve 
almost merged with the purple one). Locations of peaks as given 
by the approximative maximal LOS velocities (Sect. [5] point [5} 
given by Eq. ( 17 1 are plotted with orange lines. Profiles do not 



include stars of the host galaxy that are not part of the shell sys- 
tem and are normalized, so that the total flux equals to one. For 
parameters of the host galaxy potential see Sect. 6. 1 



2. Using the approximative LOSVD (purple curves). For the 
given shell at the chosen projected radius, Eq. ( 15 1 is a func- 
tion of only two parameters, the circular velocity v c at the 
shell edge radius and the current shell velocity v s . Assuming 
a behavior of shell brightness as a function of the shell ra- 
dius, Eq. (|T5j allows us to plot the whole LOSVD (Sect.[42i. 
However, computing the LOSVD and the peaks' positions 
requires a numerical approach in this framework. 

3. Using the approximative maximal LOS velocities (orange 
curves). Eq. ( fT7] i supplies the positions of the peaks directly. 
It differs from the previous approximation in the assump- 
tion about the galactocentric radius r,, max , from which comes 
the contribution to the LOSVD at the maximal speed. The 
assumption is that r vm . dx is given by Eq. ( 16 1, which was de- 



rived by MK98 for a stationary shell. This equation is ac- 
tually only very approximate, but allows us to analytically 
invert Eq. ( 17 1 to obtain formulae for the calculation of v c 



and v s from the measured peak positions in the spectrum of 
the shell galaxy near the shell edge (Eqs. ([T8j and (fT9]l). For 
a moving shell, we could not derive a more accurate formula 
for r vmax that would be simple enough to make the calcula- 
tion of v c and v s feasible. 

Fig. [6] shows a comparison of the radii that contribute to 
the LOSVD at the maximal velocities according to all three ap- 
proaches. For the first two methods, the radius corresponding to 
the inner maxima of the LOSVD (which are the maxima cre- 
ated by the inward stars) is lower than that for the outer maxima, 
whereas Eq. ( JT6| assumes the same r vmax for both inward and 
outward stars. 

Fig. [7] shows locations of the LOSVD peaks for the second 
shell at the radius of 120 kpc near the shell edge radius. The pur- 
ple curve is calculated using the approximative LOSVD (Sect. [5] 
point [2]i given by Eq. ( [15) , into which we inserted the velocity 
of the second shell according to the model of radial oscillations 
and the circular velocity in the potential of the host galaxy (see 



Sect. 6.1 for parameters of the potential). The purple curve does 
not differ significantly from the light blue curve calculated in the 
model of radial oscillations (Sect. 2.3 1. The more important de- 
viations in the orange curve of the approximative maximal LOS 
velocities (Sect. [5] point[3]l given by Eq. ( 17 1, are caused by us- 
ing Eq. ( 16 1 for r vmax . With this assumption, approximative max- 
imal LOS velocities (the orange curve) predict that around the 
zone between the current turning point and the shell edge, the in- 
ner peaks change signs. This means that for the part of the galaxy 
closer to the observer, both inner and outer peaks will fall into 
negative values of the LOS velocity and vice versa. However, 
from the model of the radial oscillations we know that the signal 
from the inner peak in a given (near or far) part of the galaxy is 
always zero or has the opposite sign to that of the outer peak. 

The model of the radial oscillations and the approxima- 
tive LOSVD given by Eq. ( JT5] l were also used to construct the 
LOSVD for the second shell located at 120 kpc, at the projected 
radius of 108 kpc in Fig. [8] The graph also shows the locations 
of the peaks using the approximative maximal LOS velocities 
given by Eq. ( 17 1. 



6. Test-particle simulation of the merger 

We performed a simplified simulation of formation of shells in 
a radial galactic minor merge. Both merging galaxies are repre- 
sented by smooth potential. Millions of test particles were gen- 
erated so that they follow the distribution function of the canni- 
balized galaxy at the beginning of the simulation. The particles 
then move according to the sum of the gravitational potentials 
of both galaxies. When the centers of the galaxies pass through 
each other, the potential of the cannibalized galaxy is suddenly 
switched off and the particles continue to move only in the fixed 
potential of the host galaxy. We use the simulation to demon- 
strate the validity of our methods of recovering the parameters 
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Fig. 9. Snapshots from our test-particle simulation of the radial 
minor merger, leading to the formation of shells. Each panel cov- 
ers 300x300 kpc and is centered on the host galaxy. Only the 
surface density of particles originally belonging to the satellite 
galaxy is displayed. The density scale varies between frames, so 
that the respective range of densities is optimally covered. 

of the host galaxy potential by measuring^] the positions of the 
peaks in the spectral lines. In all cases, we look at the galaxy 
from the view perpendicular to the axis of collision, so that the 
cannibalized galaxy originally flew in from the right. 

6.1. Parameters of the simulation 

The potential of the host galaxy is modeled as a double Plummer 
sphere with respective masses M, = 2 x 10" M0 and Mdm = 
1.2x10 Mq, and Plummer radii b* = 5 kpc and £>dm = 100 kpc 
for the luminous component and the dark halo, respectively. This 
model has properties consistent with observed massive early- 
type (and even s hell) galaxies (|Auger et a~ 2010 Nagino & 



Matsushita 2009 Fukazawa et al. 



cannibalized galaxy is chosen to 



|2006) . The potential of the 



*>e a single Plummer sphere 



with the total mass M — 2 x 10 M© and Plummer radius 
b t = 2 kpc. 

The cannibalized galaxy is released from rest at a distance of 
100 kpc from the center of the host galaxy. When it reaches the 
center of the host galaxy in 306.4 Myr, its potential is switched 
off and its particles begin to oscillate freely in the host galaxy. 
The shells start appearing visibly from about 50 kpc of galac- 
tocentric distance and disappear at around 200 kpc, as there are 
very few particles with apocenters outside these radii (Fig. [9}. 

6.2. Comparison of the simulation with models 

In the simulations, some of the assumptions that we used earlier 
(Sect. [2} are not fulfilled. First, the particles do not move radially, 
but on more general trajectories, which are, even in the case of 




41.111 ~ 




-100 - 50 50 
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Fig. 10. Simulated shell structure 2.2 Gyr after the decay of the 
cannibalized galaxy. Only the particles originally belonging to 
the cannibalized galaxy are taken into account. Top: surface 
density map; middle: the LOSVD density map of particles in 
the +1 kpc band around the collision axis; bottom: histogram of 
galactocentric distances of particles. The angle between the ra- 
dial position vector of the particle and the x-axis (the collision 
axis) is less than 90° for the blue curve and less than 45° for 
the red curve. The horizontal axis corresponds to the projected 
distance X in the upper panel, to the projected radius R in the 
middle panel, and to the galactocentric distance r in the lower 
panel. 



a radial merger, nevertheless very eccentric. Second, not all the 
particles are released from the cannibalized galaxy right in the 
center of the host galaxy; when the potential is switched off, the 
particles are located in the broad surroundings of the center and 
some are even released before the decay of the galaxy. These 
effects cause a smearing of the kinematical imprint of shells, as 
the turning points are not at a sharply defined radius, but rather 
in some interval of radii for a given time. 

The model of radial oscillations presented in Sect. |2]predicts 
that 2.2 Gyr after the decay of the cannibalized galaxy (Fig.[T0|, 
five outermost shells should lie at the radii of 257.3, -157.8, 
105.1, -70.5, and 48.8 kpc. The negative radii refer to the shell 
being on the opposite side of the host galaxy with respect to 
the direction from which the cannibalized galaxy flew in. These 
radii agree well with the radii of the shells measured in the sim- 
ulation 2.2 Gyr after the decay of the cannibalized galaxy. In the 
simulation, the first shell at 257.4 kpc is composed of only a few 
particles, and therefore we will not consider it (its parameters are 
listed in Table[2]for completeness). Thus, the outermost relevant 
shell in the system lies at -157.8 kpc and has a serial number 
n = 2. Also, the shell at 48.8 kpc suffers from lack of particles, 
but we will include it nevertheless. 
Fig. 



2 Here and in the rest of this section, the data measured are the output 
of our simulation. 



1 1 shows the comparison between the LOSVD in the 
simulation, the peaks of the LOSVD computed in the model 
of radial oscillations (light blue curves), and the approxima- 
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Table 2. Parameters of the shells in a simulation 2.2 Gyr after 
the decay of the cannibalized galaxy. 



r s 


n 


r TP.modcl 




^s.model 


^c.model 


kpc 




kpc 


km/s 


km/s 


km/s 


48.8 


5 


48.5 


38.7+2.1 


38.7 


326 


-70.6 


4 


-69.9 


59.8±1.6 


54.3 


390 


105.0 


3 


103.9 


68.1+1.9 


63.5 


441 


-157.8 


2 


-155.7 


74.3+1.2 


72.4 


450 


257.4 


1 


251.0 


97.5+1.4 


95.7 


406 



Notes. The values of r TPmode i and v s>mo d e i are calculated for the shell 
position r s and its corresponding serial number n according to the model 
of radial oscillations (Sect. [2J. The shell velocity v s sim is derived from 
20 positions between the times 2.49-2.51 Gyr for each shell. The value 
v c ,modei corresponds to the circular velocity at the shell edge radius r s for 
the chosen potential of the host galaxy (Sect.|6~T}, 



E 200 




6.3. Recovering the potential from the simulated data 

We used a snapshot from our simulation, which 2.2 Gyr after the 
decay of the cannibalized galaxy, as a source of the simulated 
data and tried to reconstruct the parameters of the potential of the 
host galaxy from the locations of the LOSVD peaks measured 
from the simulated data. 

For a given host galaxy, the signal-to-noise (S /N) ratio in the 
simulated data is a function of the number of simulated particles, 
the age of the shell system, the distribution function of the canni- 
balized galaxy, and the impact velocity. For a given radius in the 
simulated data, we can obtain arbitrarily good or bad S/N ratios 
by tuning these parameters. Thus, we adopted the universal cri- 
teria: 1) the LOSVD of each shell is observed down to 0.9 times 
its radius; 2) we measured the positions of the LOSVD peaks in 
different locations within the shell, sampled by 1 kpc steps. We 
do not estimate the errors, since the real data will be dominated 
by other sources. We quote only the mean square deviation and 
the standard error of the linear regression. 

These criteria give us between 7 and 15 measurements for 
a shell. Each measurement contains two values: the positions of 
the outer and inner peaks, vi os . max+ and vi os .max-, respectively, for 
each projected radius R (see green crosses in Fig. 13 1. 

First we used the approximative maximal LOS velocities 
given by Eqs. (|T~8]> and ( 19 1 for a direct calculation of the cir- 



-100 -50 



cular velocity v c eq {T8j at the shell edge radius r s and the cur- 
rent shell velocity v s , eq {T9j. These equations are the inverse of 
Eq. ( 17 1, which corresponds to the model shown in orange lines 
in pictures throughout the text (Sect. [5 point B). Mean values 
from all the measurements for each shell are shown in Tables [3] 
andgl 

We obtain a better agreement with the circular velocity of our 
host galaxy potential when using the slope of the LOSVD inten- 
sity maxima given by Eq. (21 1, where we fit the linear function 



Fig. 11. LOSVD map of the simulated shell structure 2.2 Gyr 
after the decay of the cannibalized galaxy (middle panel in 
Fig. [TO}. Light blue curves show locations of the maxima accord- 
ing to the model of radial oscillations (Sect. 2.3 1 for shell radius 
r s , corresponding serial number n, and the known potential of the 
host galaxy (Sect. |6.1). Orange curves are derived from the ap- 



of the measured distance between the outer and the inner peak 
on the projected radius (v CiS io P e in Tableland in Fig. 14 1. 

From the approximative maximal LOS velocities (Sect. [5] 
point [3 i given by Eq. ( [17} , we can derive a hybrid relation be- 
tween the positions of the LOSVD peaks, the circular velocity at 
the shell edge radius v c , and the shell velocity: 



proximative maximal LOS velocities (Sect. [5] point [3} given by vi = v l c (l-R/r s0 ) + 



^los,max+^los,max- 



17 1 for r s , v Simo dei, and v c>mo dei- Parameters of the shells are 
shown in Table [2] Black lines mark the location at 0.9r s for each 

'The 



4(fl/r s0 ) 2 (l+tf/r s0 r 2 -l 



(22) 



12 



shell. The LOSVD for these locations are shown in Fig. 
map includes only stars originally belonging to the cannibalized 
galaxy. 



tive maximal LOS velocities - Eq. ( [17} (orange curves). To eval- 
uate the approximative maximal LOS velocities, we obtained 
the shell velocity v s ,modei from the model of radial oscillations 
(Sect. [2} for the respective serial number n of the shell and circu- 
lar velocity v Cjmo dei at the shell edge radius, using our knowledge 
of the potential of the host galaxy (see Sect. 6. 1 for parameters 
of the potential). The values of all the respective shell quantities 
are listed in Table [2] 

Fig. 



11 also shows the locations that correspond to the radii 



of 0.9r s for each individual shell (black lines). The LOSVD for 



these locations is shown in Fig. 12 The positions of simulated 
LOSVD peaks largely agree with three theoretical approaches 
described in Sect. [5] 



We substitute the values of v C;S i op e derived from the measure- 
ments (that we know better describe the real circular velocity 
of host galaxy) into this relation, thus obtaining the improved 
measured shell velocity v s ,eq|22j-siope (Tableland Fig. 15 i. 

In the zone between the current turning points and the shell 
edge, the inner peaks coalesce and gradually disappear (Fig.|4}. 
The simulated data do not show a disappearance of the inner 
peaks as abrupt and clear as the theoretical LOSVD profiles pre- 
dict, so that in this zone, we can usually measure one inner peak 
at km/s. The information from these measurements is degener- 
ate, and thus we defined a subsample of simulated measurements 
with all four clear peaks in the LOSVD (in the columns labeled 
SS inTables[3]and[4}. 

The spread of the values derived using the approximative 
maximal LOS velocities given by Eqs. (18i and ([19} is signif- 
icantly lower for the subsample (v^ q rfgj and v^ q rrjj) due to the 
exclusion of areas where these equations do not nold well. On 
the contrary, the slope of the linear regression in Eq. pT} using 
the slope of the LOSVD intensity maxima gives a worse result 
(with a larger error) for the subsample . 
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Fig. 12. LOS VDs of four shells at projected radii 0.9r s (indicated as the title of each plot) 2.2 Gyr after the decay of the cannibalized 
galaxy (parameters of the shells are shown in Table [2}. The simulated data are shown in green, the LOSVDs according to the 
approximative LOSVD (Sect. [5] point |2|i given by Eq. ( 15 1 in purple, and LOSVDs according to the model of radial oscillations 
(Sect. 2.3 i in light blue. The graph also shows the locations of the peaks using the approximative maximal LOS velocities (Sect. 5] 
point|3| given by Eq. ( fTT) by orange lines. Profiles do not include stars of the host galaxy, which are not part of the shell system. The 
theoretical profiles are scaled so that the intensity of their highest peak approximately agrees with the highest peak of the simulated 
data. Intensity is given in relative units, so maxima of the profiles have values of about 0.9. 



Table 3. Circular velocity at the shell edge radius r s derived from the measurement of the simulated data 2.2 Gyr after the decay of 
the cannibalized galaxy. 



r s 

kpc 


"c^nodel 

km/s 


A'dala 


iwSS 
"data 


v c,eqp8] 

km/s 


ss 

V c,eqfi8] 

km/s 


^c.slope 

km/s 


ss 

c.slope 

km/s 


Vc.flt 

km/s 


V c ,slope(MK98) 

km/s 


48.8 


326 


5 


4 


346±130 


340±94 


322+19 


314+32 


318+51 


449+26 


-70.6 


390 


7 


5 


394±85 


390±53 


391 ±5 


392+11 


368+60 


570+23 


105.0 


441 


11 


8 


478+144 


452±64 


440±5 


447+7 


427+28 


632+9 


-157.8 


450 


15 


10 


497±236 


472±79 


462±8 


484+14 


460+32 


671+11 



Notes. r s and v c , m0 dei have the same meaning as in Table[2 A^iata^ number of measurements for each shell; v Mt| rnn: the mean of values derived from 
the approximative maximal LOS velocities given by Eq. ■ T8l with its mean square deviation; v c slope : a value derived from linear regression using 
the slope of the LOSVD intensity maxima given by Eq. i and its standard error (see also Fig.|14[l; v^,: a value derived by fitting a pair of v c 
and v s in the approximative LOSVD given by Eq. <| 1 5[ > (Sect. B] point|2]and Fig.|13|; v c . s i ope ( MK a8) '■ the mean of values derived from the slope of 
the LOSVD intensity maxima given by Eq. \2\\ with its standard error (see also Fig.fTjV In the equation, however, Avi os is substituted with the 
distance between the two outer peaks of the LOSVD intensity maxima in order to mimic the measurement as originally proposed by MK98|for 
double-peaked profile. The quantities with the superscript SS correspond to the subsample, where only measurements with two discernible inner 
peaks in the LOSVD are used. 



Table 4. Velocity of the shell at the radius r s derived from the measurement of the simulated data 2.2 Gyr after the decay of the 
cannibalized galaxy. 



r 8 

kpc 


^s.model 

km/s 


V d a(a 




km/s 


v s.eqfi9l 

km/s 


,,ss 

V S ,eqp] 

km/s 


v s,eq|22)-slope 

km/s 


v ss 

s,eq|22j— slope 

km/s 


Vs.flt 

km/s 


48.8 


38.7 


5 


4 


38.7+2.1 


50.7+2.3 


51.7+1.1 


44.2+6.5 


44.9+6.3 


53+16 


-70.6 


54.3 


7 


5 


59.8+1.6 


60.8+9.8 


65.6+2.0 


60.7+10.8 


66.0+2.9 


66+19 


105.0 


63.5 


11 


8 


68.1+1.9 


74.8+4.6 


76.5+1.4 


68.0+8.9 


71.3+2.5 


79+9 


-157.8 


72.4 


15 


10 


74.3+1.2 


84.4+5.4 


86.7+2.0 


78.7+10.5 


82.±3.5 


85+14 



Notes. r s , v s model , and v s . s j m have the same meaning as in Table [2| N^ : number of measurements for each shell; v s cq jT9]: the mean of values derived 
from the approximative maximal LOS velocities given by Eq!~|l9[l with its mean square deviation; v s>cq [22]-siope: the mean of values derived from 
the hybrid relation given by Eq. l |22[ ) with its mean square deviation (see also Fig.|15}; v s fl t : a value derived by fitting a pair of v c and v s in the 
approximative LOSVD given by Eq. \\5) (Sect.[5] point|2]and Fig.|13|>. The quantities with the superscript SS correspond to the subsample, where 
only measurements with two discernible inner peaks in the LOSVD are used. 



The third option to derive the circular velocity v c at the shell 
edge radius r s and shell velocity v s from the simulated data is to 
use the approximative LOSVD given by Eq. ( 15 1, which corre- 
sponds to the model shown in purple lines in pictures throughout 
the text (Sect. [5] point[2]i. However, this requires a numerical so- 
lution of the equation for a given pair of v c and v s . We minimized 
the sum of sums of squared differences between vi os max± (v c , v s ) 
as given by the approximative LOSVD and the simulated data 



to obtain best fitted values v Cj fi t and v Sj fi t (see Tables [3] and |4] 
for the results). Errors were estimated using the ordinary least 
squared minimization as if the functions vi s,max+(v c ,fit, fs.fit) and 
vios,max-(v c ,fit, v s,fit) were fitted separately; quoted is the larger of 
the two errors. The LOSVD intensity maxima resulting from this 



procedure are plotted in Fig. 13 together with the fitted data and 



the maxima given by the model of radial oscillations (Sect.|23 
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Fig. 13. Fits for circular velocity v c and shell velocity v s using the approximative LOSVD (Sect, i] point [ij given by Eq. ( 15 1 for 
four shells (r s indicated in bottom right corner of each plot) in the simulation 2.2 Gyr after the decay of the cannibalized galaxy. 
The best fit is the purple curve, and its parameters are shown in Tables [3] and [4] in the columns labeled v Cj fit and v S( fl t . The green 
crosses mark the measured maxima in the LOSVD, and the light blue curves show the locations of the theoretical maxima derived 
from the host galaxy potential according to the model of radial oscillations (Sect. |2~3| l. Note that the values of v c and v s used in the 
approximative LOSVD for the purple line were obtained by fitting the parameters to the simulated data, whereas in Figs. [7] [8} and 
[12] the values are known from the model of the host galaxy potential. 



For the sake of comparison with the method of MK98 we 
calculated the circular velocity v CiS i ope (MK98) at the shell edge ra- 
dius r s using the slope of the LOSVD intensity maxima given 
by Eq. (21 1. To mimic the measurement of the circular veloc- 



ity according to the Eq. (7) in MK98 which was derived for the 
double-peaked profile, we assume Avi os is the distance between 
the two outer peaks of the LOSVD intensity maxima. In Table[3] 



and Fig. 14 



we can easily see that the values v c s iope(MK98) differ 
from the~actual circular velocity of the host galaxy v c m0 dei by a 
factor of 1.3-1.5. 



6.4. Instrumental LOSVD 



When observed, the LOSVD is always influenced by instrumen- 
tal dispersion, which naturally smoothes features of the spectral 
profile. In Fig. 16 we show the LOS VDs from the simulated data 



smoothed with different Gaussians representing the instrumental 
profiles having the full width at half maximum (FWHM) of 10, 
30, and 60km/s. It is obvious that relatively high spectral reso- 
lution is necessary for observing an imprint of shell peaks in line 
profiles. 



700 



600 




host galaxy model - 

v c,slope(MK98) I- 



90 
r [kpc] 



120 



150 



180 



Fig. 14. Circular velocity of the model and values derived from 
the simulated data: v c , mo dei of the host galaxy model is shown 
by the black line; blue and red points show values of circular 
velocity as they result from the analysis of the simulated LOSVD 
(see Sect. |6.2| and Table|3]for the numbers). 
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Fig. 16. Line profiles of four shells at projected radii 0.9r s (indicated as the title of each plot, same as in Fig. 12 1 2.2 Gyr after the 
decay of the cannibalized galaxy: gray lines show the LOSVDs for the host galaxy at a given radius (except for the radius of 44 kpc 
the signal of the host galaxy is negligible comparing to the signal from the cannibalized galaxy); green lines show the total LOSVDs 
from the host and the cannibalized galaxy together; red, blue, and yellow lines show convolutions of the total simulated data with 
different Gaussians representing the instrumental profiles having the FWHM 10, 30, and 60 km/s, respectively. Scaling is relative, 
similar as in Fig. 12 
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v TP (n = 1), host galaxy model 
[3/(2n+1)] v s , eq(2 2)-slope 
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90 
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Fig. 15. Comparison of velocity of the shell as a function of ra- 
dius from the model and the simulated data. Velocity for the first 
shell (n = 1) in the host galaxy model is shown by the black line. 
Red crosses show V Sj eqj22l-siope (Table 4) as they result from the 
analysis of the simulateaLOSVD. Values are corrected for shell 
number n by the factor 3/(2n + 1), so they correspond to velocity 
of the first shell (e.g., Eq. (3]l). 



7. Discussion 

We developed a new method to measure the potential of shell 
galaxies from kinematical data, extending the work of MK98 
The method splits into three different analytical and semi- 
analytical approaches for obtaining the circular velocity in the 
host galaxy, v c , and the current shell phase velocity, v s : 

1. The approximative LOSVD: using Eq. (JT3J and an assump- 
tion of the behavior of the shell brightness as a function of 
the shell radius (Sect. |4~2| l. 

2. The approximative maximal LOS velocities: Eqs. (T8} and 
(I9j (Sectgg). 

3. Using the slope of the LOSVD intensity maxima in the 
R x vjos diagram in Eq. pT) ) (Sect. |4.4[ ). 



In Sect. [5] the first two approaches are compared to the model 
of radial oscillations (numerical integration of radial trajecto- 
ries of stars in the host galaxy potential, Sect. [2]). All three ap- 
proaches are then applied to data for the four shells obtained 
from a test-particle simulation and compared to the theoretical 
values (Sect. [oT2| i. 

Approach[l]requires a numerical solution to Eq. (T5) and the 
search for a pair of v c and v s , which matches the (simulated) data 
best. Although this approach is not limited by any assumptions 



about the radius of the maximal LOS velocity (Sect. 4.3 I, it does 
not give a better estimate of v c and v s for our simulated shell 
galaxy than the other two methods. The deviation from the real 
value of v c is between 2 % and 6 %. 

Using the approximative maximal LOS velocities ap- 
proach [2] results in simple analytical relations and is the only 
one that can in principle be used for an LOSVD measured at 
only one projected radius. Nevertheless, when measuring in the 
zone between the radius of the current turning points and the 
shell radius, we can expect very bad estimates of v c and v s . The 
mean value from more measurements of the LOSVD peaks for 
each shell of our simulated shell galaxy has similar accuracy to 
those of approach[TJ provided that we include only the measure- 
ments outside the zone between the radius of the current turning 
points and the shell radius. 

The best method for deriving the circular velocity in the po- 
tential of the host galaxy seems to be to use the slope of the 
LOSVD intensity maxima, with a typical deviation in the order 
of units of km/ s when fitting a linear function over all the mea- 
sured positions of the LOSVD peaks for each shell. This circular 
velocity is then used in the hybrid relation, Eq. (22]), to obtain the 
best estimate of the shell velocity. 

All the approaches, however, derive a shell velocity system- 
atically larger than the prediction of the model of radial oscil- 
lations v s mo( iei and the value derived from positions between the 
times 2.49-2.51 Gyr in the simulation v s s ; m (Table[3]l. This is be- 
cause the simulated LOSVD peaks lie too far out (for the outer 
peaks) or too far in (for the inner peaks) when compared to the 
model of radial oscillations. That can be caused by nonradial 
trajectories of the stars of the cannibalized galaxy or by poor 
definition of the shell radius in the simulation. 
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Nevertheless, the shell velocity depends, even in the simpli- 
fied model of an instant decay of the cannibalized galaxy in a 
spherically symmetric host galaxy (Sec. |2j, on the serial number 
of the shell n and on the whole potential from the center of the 
galaxy up to the shell radius (see Eq. Q). A comparison of its 
measured velocity to theoretical predictions is possible only for a 
given model of the potential of the host galaxy and the presumed 
serial number of the observed shells. In such a case, however, it 
can be used to exclude some parameters or models of the poten- 
tial that would otherwise fit the observed circular velocity. 

The first shell has a serial number equal to one. A higher 
serial number means a younger shell. On the same radius, the 
velocity of each shell is always smaller than that of the previous 
one. In practice, it is difficult to establish whether the outermost 
observed shell was the first one created, or whether the first shell 
(or even the first couple of shells) were already unobservable. 
Here, we can use the potential derived from our method or a 
completely different one in a reverse way: to determine the ve- 
locity of the first shell on the given radius and to compare it to 
the velocity derived from the positions of the LOSVD peaks. 
Knowing the serial number of the outermost shell and its posi- 
tion allows us then to determine the time from the merger and the 
impact direction of the cannibalized galaxy. Moreover, the mea- 
surement of shell velocities reveals the presence of shells from 
different generations (Bartos kova et al.|201 1 1. 

Our method for measuring the potential of shell galaxies has 
several limitations. Theoretical analyses were conducted over 
spherically symmetric shells, while the test-particle simulation 
was run for a strictly radial merger and analyzed in a projection 
plane parallel to the axis of the merger. In addition, both ana- 
lytical analysis and simulations assume spherical symmetry of 
the potential of the host galaxy. In reality, the regular shell sys- 
tems with more shells in one galaxy are more often connected 
to galaxies with significant elfipticity (Dupraz & Combes 1986). 
Moreover, in cosmological simulations with cold dark matter, 
halos of galaxies are described as triaxial ellipsoids (e.g., |Jing| 



& Suto 2002; Bailin & Steinmetz 2005; Allgood et al. 2006). 

ii " i ii 1 1 - j i i 



However, the effect of the ellipticity of the isophotes of the host 
galaxy on the shell kinematics need not be dramatic, as the shells 
have the tendency to follow equipotentials that are in general 
less elliptical than the isophotes. Dupraz & Combes ([1986 ) con- 
cluded that while the ellipticity of observed shells is generally 
low, it is neatly correlated to the eccentricity of the host galaxy. 



Prieur ( 1988 1 pointed out that the shells in NGC 3923 are much 



rounder than the underlying galaxy and have an ellipticity that 
is similar to the inferred equipotential surfaces. This idea was 
originally put forward by Dupraz & Combes ( 1986 1, who found 
such a relationship for their merger simulations. Our method is 
in principle applicable even to shells spread around the galactic 
center, which are usually connected to rounder elliptical galaxies 
if they were created in a close-to-radial merger. Nevertheless, the 
combination of the effects of the projection plane, merger axis, 
and ellipticity of the host galaxy can modify our results and re- 
quire further analyses. 

Because the kinematics of the stars that left the cannibalized 
galaxy is in the first approximation a test-particle problem, they 
should not be much affected by self-gravity of the cannibalized 
galaxy and the dynamical friction that this galaxy undergoes dur- 
ing the merger, both of which have been neglected in this work. 

Another complication is that the spectral resolution required 



to distinguish all four peaks is probably quite high (Sect. 6.4 



and Fig.[16[i and the shell contrast is usually small. Nevertheless, 
there is the possibility to measure shell kinematics using the LOS 



velocities of individual globular clusters, planetary nebulae, and, 
in the Local group of galaxies, even of individual stars. 

8. Conclusions 

Kinematics of regular shells produced during nearly radial mi- 
nor mergers of galaxies can be used to constrain their gravita- 
tional force field and thus the dark matter distribution. Merrifield 
|& Kuijken] ( |1998[ ) showed that the LOSVD measured near the 
edges of a shell has a double-peaked shape, and found a relation 
between the values of the two LOS velocity peaks and the circu- 
lar velocity. Their approximation is limited to stationary shells. 

We have extended their theoretical analysis to traveling 
shells. We find that in two-component giant galaxies with real- 
istically massive dark matter halos, shell propagation velocity is 
significantly higher, typically 30-150 km/s, compared to values 
quoted in the theoretical studies in the literature. We show that 
such large speeds have considerable impact on the LOS kine- 
matics of shells. We demonstrate that each peak of the double- 
peaked profile is split into two, producing a quadruple-peaked 
LOSVD. We derive a new approximation, relating the circular 
velocity of the host galaxy potential at the shell edge radius, as 
well as the current phase velocity of the shell, to the positions of 
the four peaks. 

In galaxies with multiple shells, we can use circular veloc- 
ities measured by these methods to determine the potential of 
the host galaxy over a large span in radii, whereas the measured 
shell phase velocity carries information on the age of the shell 
system, and the arrival direction of the cannibalized galaxy. The 
potential observation of multigeneration shell systems contains 
additional limits on the shape of the potential of the host galaxy. 

Acknowledgements. We acknowledge support from the following sources: grant 
No. 205/08/H005 by Czech Science Foundation (IE, LJ, MK, KB, and TS); grant 
MUNI/A/0968/2009 by Masaryk University in Brno (LJ, KB, and TS); research 
plan AV0Z 10030501 by Academy of Sciences of the Czech Republic (IE, BJ, 
MK, MB, KB, and IS); grant LC06014- Center for Theoretical Astrophysics - 
by Czech Ministry of Education (BJ and IS); and the project SVV261301 
by Charles University in Prague (IE, MK, and MB). This work has been 
done with the support for a long-term development of the research institution 
RV067985815 (IE, BJ, MB, and KB). LJ acknowledges the support by the 2- 
years ESO PhD studentship, held in ESO, Santiago. 



References 

Allgood, B., Flores, R. A., Primack, J. R., et al. 2006, MNRAS, 367, 1781 
Arp, H. 1966, ApJS, 14, 1 

Auger, M. W., Treu, T., Bolton, A. S., et al. 2010, ApJ, 724, 511 
Bailin, J. & Steinmetz, M. 2005, ApJ, 627, 647 
Balcells, M. & Sancisi, R. 1996, AJ, 111, 1053 

Bartoskova, K, Jungwiert, B., Ebrova, I., Jflkova, L., & Kffzek, M. 2011, 
Simulations of Shell Galaxies with GADGET-2: Multi-Generation Shell 
Systems, ed. I. Ferreras & A. Pasquali, 195 
Canalizo, G., Bennert, N., Jungwiert, B., et al. 2007, ApJ, 669, 801 
Charmandaris, V., Combes, E, & van der Hulst, J. M. 2000, A&A, 356, LI 
Churazov, E., Forman, W., Vikhlinin, A., et al. 2008, MNRAS, 388, 1062 
Churazov, E., Tremaine, S., Forman, W., et al. 2010, MNRAS, 404, 1165 
Coccato, L., Gerhard, O., Arnaboldi, M., et al. 2009, MNRAS, 394, 1249 
Colbert, J. W., Mulchaey, J. S., & Zabludoff, A. I. 2001, AJ, 121, 808 
Coleman, M., Da Costa, G. S., Bland-Hawthorn, J., et al. 2004, AJ, 127, 832 
Das, P., Gerhard, O., Churazov, E., & Zhuravleva, I. 2010, MNRAS, 409, 1362 
de Lorenzi, E, Gerhard, O., Coccato, L., et al. 2009, MNRAS, 395, 76 
Deason, A. J., Belokurov, V., Evans, N. W., & McCarthy, I. G. 2012, ApJ, 748, 
2 

Due, P.-A., Cuillandre, J.-C, Serra, P., et al. 201 1, MNRAS, 417, 863 
Dupraz, C. & Combes, F. 1986, A&A, 166, 53 
Dupraz, C. & Combes, F. 1987, A&A, 185, LI 

Ebrova, I., Jflkova, L., Jungwiert, B., et al. 2011, Quadruple-Peaked Line-of- 
Sight Velocity Distributions in Shell Galaxies, ed. I. Ferreras & A. Pasquali, 
225 



14 



I. Ebrova et al.: Quadruple-peaked spectral line profiles as a tool to constrain gravitational potential of shell galaxies 



Ebrova, I., Jungwiert, B., Canalizo, G., Bennert, N., & Jflkova, L. 2010, in 
Astronomical Society of the Pacific Conference Series, Vol. 423, Galaxy 
Wars: Stellar Populations and Star Formation in Interacting Galaxies, ed. 

B. Smith, J. Higdon, S. Higdon, & N. Bastian, 236 

Fardal, M. A., Babul, A., Guhathakurta, P., Gilbert, K. M., & Dodge, C 2008, 
ApJ, 682, L33 

Fardal, M. A., Guhathakurta, P., Babul, A., & McConnachie, A. W. 2007, 
MNRAS, 380, 15 

Fardal, M. A., Guhathakurta, P., Gilbert, K. M., et al. 2012, MNRAS, 3140 
Fukazawa, Y., Botoya-Nonesa, J. G, Pu, J., Ohto, A., & Kawano, N. 2006, ApJ, 
636, 698 

Gavazzi, R., Treu, T., Rhodes, J. D., et al. 2007, ApJ, 667, 176 
Heisler, J. & White, S. D. M. 1990, MNRAS, 243, 199 

Helmi, A., Navarro, J. F., Meza, A., Steinmetz, M., & Eke, V. R. 2003, ApJ, 592, 
L25 

Hemquist, L. & Quinn, P. J. 1987a, ApJ, 312, 17 

Hernquist, L. & Quinn, P. J. 1987b, ApJ, 312, 1 

Hemquist, L. & Quinn, P. J. 1988, ApJ, 331, 682 

Hernquist, L. & Spergel, D. N. 1992, ApJ, 399, LI 17 

Horellou, C, Black, J. H., van Gorkom, J. H., et al. 2001, A&A, 376, 837 

James, R. A. & Wilkinson, A. 1987, in IAU Symposium, Vol. 127, Structure and 

Dynamics of Elliptical Galaxies, ed. P. T. de Zeeuw, 471 — t 
Jflkova, L., Jungwiert, B., Kfizek, M., et al. 2010, in Astronomical Society of 

the Pacific Conference Series, Vol. 423, Galaxy Wars: Stellar Populations and 

Star Formation in Interacting Galaxies, ed. B. Smith, J. Higdon, S. Higdon, 

& N. Bastian, 243 
Jing, Y. P. & Suto, Y. 2002, ApJ, 574, 538 
Kaviraj, S. 2010, MNRAS, 406, 382 
Kim, T., Sheth, K., Hinz, J. L., et al. 2012, ApJ, 753, 43 
Krajnovic, D., Emsellem, E., Cappellari, M., et al. 201 1, MNRAS, 414, 2923 
Malin, D. F. & Carter, D. 1983, ApJ, 274, 534 

Mandelbaum, R., Seljak, U., & Hirata, C. M. 2008, J. Cosmology Astropart. 
Phys., 8, 6 

Marino, A., Iodice, E., Tantalo, R., et al. 2009, A&A, 508, 1235 
Merrifield, M. R. & Kuijken, K. 1998, MNRAS, 297, 1292 
Miskolczi, A., Bomans, D. J., & Dettmar, R.-J. 2011, A&A, 536, A66 
Nagino, R. & Matsushita, K. 2009, A&A, 501, 157 

Nierenberg, A. M., Auger, M. W., Treu, T., Marshall, P. J., & Fassnacht, C. D. 
2011, ApJ, 731,44 

Norris, M. A., Gebhardt, K., Sharpies, R. M., et al. 2012, MNRAS, 421, 1485 
Nulsen, P. E. J. 1989, ApJ, 346, 690 

Petric, A., Schiminovich, D., van Gorkom, J., van der Hulst, J. M., & Weil, M. 

1997, in Bulletin of the American Astronomical Society, Vol. 29, Bulletin of 

the American Astronomical Society, 1344 — t 
Prieur, J.-L. 1988, ApJ, 326, 596 
Quinn, P. J. 1984, ApJ, 279, 596 

Rampazzo, R., Marino, A., Tantalo, R., et al. 2007, MNRAS, 381, 245 
Rampazzo, R., Plana, H., Longhetti, M., et al. 2003, MNRAS, 343, 819 
Romanowsky, A. J., Strader, J., Brodie, J. P., et al. 2012, ApJ, 748, 29 
Sanderson, R. E., Mohayaee, R., & Silk, J. 2012, MNRAS, 420, 2445 
Schiminovich, D., van Gorkom, J. H., van der Hulst, J. M., & Kasow, S. 1994, 
ApJ, 423, L101 + 

Schiminovich, D., van Gorkom, J. H., van der Hulst, J. M., & Malin, D. F. 1995, 
ApJ, 444, L77 

Schweizer, F. 1983, in IAU Symposium, Vol. 100, Internal Kinematics and 

Dynamics of Galaxies, ed. E. Athanassoula, 319-326 
Schweizer, F. & Ford, Jr., W. K. 1985, in Lecture Notes in Physics, Berlin 

Springer Verlag, Vol. 232, New Aspects of Galaxy Photometry, ed. J.-L. 

Nieto, 145-+ 

Sikkema, G, Carter, D., Peletier, R. F, et al. 2007, A&A, 467, 1011 

Stickel, M., van der Hulst, J. M., van Gorkom, J. H., Schiminovich, D., & Carilli, 

C. L. 2004, A&A, 415, 95 

Thomas, J., Saglia, R. P., Bender, R., et al. 201 1, MNRAS, 415, 545 
Trinchieri, G, Rampazzo, R., Chiosi, C, et al. 2008, A&A, 489, 85 
Weijmans, A.-M., Cappellari, M., Bacon, R., et al. 2009, MNRAS, 398, 561 
Weijmans, A.-M., Krajnovic, D., van de Ven, G, et al. 2008, MNRAS, 383, 1343 



15 



